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Описан процесс оптимизации оптических свойств плазменных теплоотражающих покрытий. Расчёт построен на использовании 
так называемого эволюционного метода, который состоит в целенаправленном переборе компонентов покрытия с последова- 
тельным приближением к глобальному экстремуму по законам биологического сообщества. Он позволяет существенно сокра- 
тить трудоёмкость процесса оптимизации. Показано, что для достижения высоких коэффициентов пропускания теплоотражаю- 
щих покрытий выгодно (для улучшения качества просветления) разделить один слой серебра на два с прослойкой просветляю- 
щих оксидных плёнок между ними. Это несколько усложняет технологию осаждения теплоотражающих покрытий, но улучшает 
их функциональные характеристики. Приведены данные по пропусканию и отражению рассмотренных покрытий. 


Введение 

Известно, что в тепловом балансе жилых зда- 
ний и промышленных объектов до 40 % потерь 
приходится на утечку тепла через прозрачные ог- 
раждения (окна, балконы, витрины и т.д.) [1]. При 


этом доля радиационной компоненты составляет 
примерно 2/3 этих потерь, и только 1/3 приходится 
на теплопроводность. 

Максимум теплового излучения жилого помеще- 
ния при 20 °С находится в инфракрасной (ИК) об- 
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ласти на длине волны 10 мкм [2]. Теплоотражающее 
покрытие (ТОП) отражая более 90 % ИК излучения, 
хорошо пропускает видимый свет. Нанесение на по- 
верхность стекла ТОП позволяет снизить потери на 
излучение в 2...3 раза по сравнению с обычным стек- 
лом, вследствие уменьшения коэффициента эмис- 
сии в 10.. .20 раз. Уменьшая потери тепла из здания в 
зимнее время, такие стекла ограничивают его неже- 
лательный приток внутрь помещения летом, благо- 
даря чему снижаются потери энергии на кондицио- 
нирование воздуха и вентиляцию. 

Раньше задача об оптимальных по материаль- 
ным и геометрическим параметрам ТОП на листо- 
вом архитектурном стекле не вызывала большого 
практического интереса из-за того, что не сущест- 
вовало технологий их равномерного нанесения на 
большие поверхности. Сегодня ситуация сущест- 
венно иная. Благодаря большим успехам в разви- 
тии плазменных методов обработки материалов эта 
проблема в общих чертах решена. Но эффектив- 
ность подобных покрытий и их потребительские 
свойства во многом зависят от того, насколько пра- 
вильно выбраны материал каждого из слоёв, их 
последовательное расположение и толщина. 

Цель настоящей работы - изложить методику 
оптимизации оптических свойств ТОП на основе 
эволюционного алгоритма, широко распростра- 
нённого в моделировании динамики развития био- 
логических сообществ. 

В данной статье однослойные образования на 
подложке мы будем называть плёнками, многос- 
лойные -покрытиями. 

Теплоотражающие покрытия как объект 

оптимизации 

Современные ТОП содержат тонкую полупроз- 
рачную плёнку металла, обладающую большим ко- 
эффициентом отражения в ИК области (обычно 
это серебро), и слои диэлектриков, чаще всего ок- 
сидов металлов, для просветления в видимой об- 
ласти спектра. 

Плёнка серебра заключается между тонкими 
(толщиной порядка 1 нм) барьерными слоями из 
коррозионно-стойкого металла или сплава. Основ- 
ная задача их состоит в том, чтобы защитить сереб- 
ряную плёнку от воздействия химически активных 
ионов кислорода, которые присутствуют в плазме 
при осаждении просветляющих (оксидных) слоёв, 
входящих в состав ТОП. Барьерные слои с одной 
стороны должны быть как можно тоньше, чтобы не 
очень сильно снижать пропускание видимого све- 
та. Но с другой стороны они должны быть доста- 
точно плотными, чтобы защитить плёнку серебра. 

Практически все технологии нанесения ТОП на 
поверхность листового стекла построены на ис- 
пользовании магнетронных распылительных сис- 
тем. Главными критериями качества стекла с ТОП 
являются его сопротивление тепловому потоку 
(чем больше отражением в ИК области, тем больше 


сопротивление тепловому потоку) и коэффициент 
пропускания света в видимой области. 

Такое стекло обычно используют в составе стек- 
лопакета покрытием внутрь. Это обеспечивает сох- 
ранность ТОП от механических воздействий и пре- 
пятствует конденсации атмосферной влаги на пок- 
рытии, а также снижает потери тепла за счет теп- 
лопроводности. 

Оптимизация ТОП заключается в варьирова- 
нии толщины просветляющих слоев до достиже- 
ния максимального пропускания в видимой облас- 
ти спектра при условии, что сопротивление тепло- 
вому потоку ТОП удовлетворяет требованиям стан- 
дарта. При этом наибольшее значение тол щ ины 
просветляющих слоев задаётся так, чтобы процесс 
интерференции не перешел в ИК область. 

Критерием качества ТОП является оценочная 
функция/ которая вычисляется по видимому диа- 
пазону спектра как сумма квадратов разностей за- 
данной и расчетной величин пропускания и/или 
отражения: 

М 

/ = ЕС 1 - (Го, - Щ.)) 2 " - Щ )) 2 ), ( 1 ) 

і=1 

где Г 0 , Го - заданные и Г, Г - расчетные величины 
коэффициентов пропускания и отражения ТОП; 
Я - длина волны падающего излучения; М - число 
точек расчета. 

Необходимо отметить, что ТОП - это многос- 
лойная структура, и величины Г и Г зависят также 
от толщины И и вида материала каждого слоя. 

Стекло с простым ТОП представляет собой ком- 
бинацию слоёв: стекло {И= 4.. .6 мм)/ТЮ 2 (20...40 нм) 
/А§(9... 13 нм)/М-Сг(~1 нм)/ТЮ 2 (20...40 нм). Изме- 
няя методом последовательных приближений зна- 
чения толщин верхнего и нижнего слоев Ті0 2 (тол- 
щина плёнки серебра берётся постоянной), можно 
добиться коэффициента пропускания видимого 
света на уровне 84... 90 %. В этом случае оценочная 
функция зависит от двух переменных и имеет один 
экстремум [3, 4]. 

По указанному выше способу нами была прове- 
дена оптимизация ТОП вида Ті0 2 /А§(11 нм)/ТЮ 2 . В 
результате измерений сопротивление тепловому по- 
току однокамерного стеклопакета, содержащего 
стекло с таким ТОП, оказалось равным 0,60 °С-Вт/м 2 , 
в то время как по требованию стандарта (СНиП 
11-3-79 "Строительная теплотехника") оно должно 
быть не ниже 0,56 °С-Вт/м 2 [1]. При этом интег- 
ральное пропускание видимого света составило 
86 % (пропускание стекла без ТОП ~90 %). 

Что надо сделать, если потребуется дополнитель- 
но увеличить сопротивление тепловому потоку в ре- 
зультате, например, ужесточения стандарта? Боль- 
шое сопротивление тепловому потоку может быть 
достигнуто простым увеличением толщины слоя 
металлической плёнки. Но это приводит к дополни- 
тельным трудностям. Например, при увеличении 
толщины слоя металла для обеспечения пропуска- 
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ния на указанном выше уровне следует использо- 
вать слои оксидов с еще более высокими показате- 
лями преломления, что очень трудно выполнить на 
практике. В этом случае один такой слой моделиру- 
ют системой чередующихся слоев оксидов с боль- 
шим и малым показателями преломления (напри- 
мер, ТіОі и АШ,), а оценочная функция начинает 
зависеть от трех и более переменных (толщин сло- 
ев). Она становится многоэкстремальной, и процесс 
оптимизации последовательным приближением ос- 
танавливается в первом локальном экстремуме 
функции, не достигая глобального экстремума. 

Возникает необходимость использования мето- 
да, который бы рассматривал все или большинство 
возможных комбинаций толщин просветляющих 
слоев. Другими словами, задача ставится так: как, 
по какому алгоритму нужно менять значения тол- 
щины просветляющих слоев в той или иной ком- 
позиции ТОП, чтобы добиться максимального эф- 
фекта просветления? 

Для её решения (при практически полном отсут- 
ствии предположений о характере исследуемой функ- 
ции) применяется так называемый эволюционный 
метод поиска, впервые предложенный Дж. Холлан- 
дом [5]. Основные положения его изложены ниже. 

Сущность эволюционного метода 

Работа эволюционного метода оптимизации 
представляет собой итерационный процесс, схема 
которого приведена на рис. 1 . Здесь в качестве объ- 
екта оптимизации выступает индивид - ТОП за- 
данной композиции на стекле. 



Рис. 1 . Схема эволюционного метода оптимизации 


Вначале формируется некое множество инди- 
видов со случайно заданными значениями толщин 
просветляющих слоев. Рассчитывается средняя 
оценочная функция: 

( 2 ) 

где Я - число индивидов. 

Затем с помощью оператора отбора из этого 
множества выбираются два индивида, у которых 
значения оценочных функции большие, чем сред- 
нее по множеству. Выбранные индивиды подверга- 
ются действию операторов кроссовера и мутации. 


Кроссовер послойно (за исключением плёнки се- 
ребра), по определенному правилу меняет значения 
толщин слоев этих индивидов. Правилом может быть 
некая формула (для десятичных цифр) или изменение 
некоторых разрядов (для двоичных цифр). Оператор 
мутации воздействует на значение тол щ ины одного из 
слоев покрытия. Он призван скачкообразно изменить 
оценочную функцию индивида в лучшую или худшую 
(в смысле оптимизации) сторону. Слишком частое 
действие оператора мутации приводит к постоянным 
скачкообразным изменениям значений оценочной 
функции. В результате средняя оценочная функция 
ведет себя непредсказуемо и процесс оптимизации 
прекращается. Поэтому действие оператора мутации 
имеет малую вероятность. 

Новое множество формируется, например, до 
тех пор, пока число индивидов в нем не будет равно 
числу индивидов в предыдущем. Далее предыдущее 
множество заменяют новым. При этом рассчитыва- 
ется новая средняя оценочная функция и т.д. 

Таким образом, средняя оценочная функция, 
увеличивая от итерации к итерации свое значение, 
задает направление изменения толщин просветля- 
ющих слоев. В результате по истечению заданного 
количества итераций остается множество индиви- 
дов с толщиной просветляющих слоев, значение 
которой близко к толщине, соответствующей гло- 
бальному экстремуму оценочной функции. 

Время счёта по алгоритму эволюционного ме- 
тода зависит от параметров и вида операторов от- 
бора, кроссовера и мутации, условий формирова- 
ния нового множества, числа индивидов в множе- 
стве, типа композиции индивида, точности расчета 
его оценочной функции что достигается увеличе- 
нием числа точек расчёта М, см. формулу (1), усло- 
вий остановки процесса оптимизации. 

Модель пропускания и отражения одно- и 

многослойной тонкой пленки 

Как видно из формулы (1), определение оценоч- 
ной функции основано на расчете Т и К. При вы- 
числении этих величин использовались выражения 
Френеля [ 6 ] для плоской изотропной пленки, тол- 
щина которой соизмерима с длиной волны падаю- 
щего света, обладающей показателем преломления 
п 2 и коэффициентом поглощения к 2 (см. рис. 2 , а, 
здесь параметры п ь к ь я 3 , к } характеризуют две по- 
лубесконечные среды по обе стороны пленки). Для 
простоты угол падения света принят нормальным. 

Обозначим (і=(2л/Х)к(п 2 +ік 2 ) и составим ряд для 
прошедшей через пленку и отраженной от пленки 
компонент падающей волны. Как видно из рис. 2, а, 
каждый ряд представляет собой бесконечную гео- 
метрическую прогрессию, членами которого явля- 
ются произведения коэффициентов Френеля г а , г 2і , 
/і 2 , /23 для границ раздела сред. Найдя сумму каждого 
ряда, получим коэффициенты отражения и пропус- 
ке ІВ 

г 12 + г 23 е г п+ г п е 

кания пленки: г 13 = 23 р/і ;г 13 = и ~ 3 ,. 2|8 ■ 

1 + П2 г 23 е 1 + '12'23 е 
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а) 



нами, константы для серебра взяты из работы [8]. 
На рис. 3 показана дисперсия оптических констант 
этих материалов в видимой области. Все они, кро- 
ме серебра, в расчётах приняты как не поглощаю- 
щие излучение свет в видимой области спектра. 
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Рис. 3. Оптические константы материалов слоев ТОП 


Рис. 2. Оптическая модель однослойной плёнки (а) и мно- 
гослойного покрытия (б) 

Величины амплитуд и фаз ср определяются сле- 
дующим образом: К=\г и \ 2 , Т = — -^1к 3 Г, 

К + 

Ііті(г п ) Іт(Г п ) 

= . \ Лёі<Рг) =77-7^- Фазы вычисляются 

К.е(л в ) К.е(Г 13 ) 

в соответствии со знаками числителя и знаменателя. 

Если на поверхности подложки имеется мно- 
гослойное покрытие (см. рис. 2, б), то расчет удоб- 
но начинать с нижней плёнки. Вычисление коэф- 
фициентов /у, и / 13 производится аналогично случаю 
однослойной пленки, находящейся между двумя 
средами. Затем поверхность с пленкой заменяется 
некой эффективной поверхностью, характеризуе- 
мой рассчитанными коэффициентами, играющи- 
ми роль френелевских коэффициентов г ъ , Г 23 . 

Этот прием повторяется до тех пор, пока не бу- 
дет добавлен последний (самый верхний) слой 
покрытия [7]. 

Особенностью стеклянной подложки является 
её очень большая по сравнению с покрытием тол- 
щина (4 мм). При попытке рассчитать для данного 
случая пропускание и отражение имеет место зна- 
чительное искажение результатов вследствие влия- 
ния многолучевой интерференции. Поэтому для 
подложки указанные параметры вычислялись без 
учета интерференции. 

Проектирование теплоотражающих покрытий. 

Результаты расчётов 

Рассмотрим ТОП на основе серебра и оксидов 
титана, алюминия и олова. Оптические константы 
оксидов и стеклянной подложки были рассчитаны 


На рис. 4 приведена расчетная зависимость интег- 
рального коэффициента пропускания видимого све- 
та серебряной пленкой от ее толщины. Видно что, по 
мере увеличения толщины плёнки от 10 до 50 нм он 
падает в 20 раз и составляет всего 3 %. Поэтому прос- 
ветление относительно плотной пленки представляет 
собой самостоятельную весьма сложную задачу. В 
данной работе она не будет рассматриваться. 



Рис. 4. Зависимость интегрального коэффициента пропус- 
кания видимого света серебряной пленки от ее тол- 
щины 

Теперь попробуем просветлить слой серебра 
толщиной 20 нм. При этом потребуем, чтобы коэф- 
фициент пропускания стекла с ТОП, содержащим 
такой слой, в видимой области спектра был бы на 
уровне стекла без каких-либо покрытий (~90 %). 

Зададим Г 0 =ЮО%иЛо = 0%. Толщина плёнки 
серебра в процессе расчета меняться не будет. Зна- 
чения толщины слоёв оксидов будем изменять в 
интервале 5... 63 нм. Условием остановки процесса 
оптимизации сделаем достижение определенного 
количества итераций, причем число индивидов в 
множестве будет оставаться фиксированным. 

Процедура отбора эволюционного метода в на- 
шем случае построена на принципе элитарности. В 
соответствии с ним лучшие индивиды без каких- 
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либо условий переходят в следующее множество. 
Таким образом, исключается потеря оптимальных 
значений толщин просветляющих слоев при пере- 
ходе от множества к множеству. 

В результате расчётов получена структура 1, схе- 
матично изображённая на рис. 5. Её коэффициен- 
ты пропускания и отражения в видимой области 
приведены на рис. 6. Высокая степень прозрачнос- 
ти (на уровне чистого стекла) достигнута лишь в 
области длин волн 500. ..600 нм. За её пределами ко- 
эффициент пропускания весьма мал. 

Дальнейшие поиски не привели к улучшению 
этих характеристик. Вначале, с увеличением коли- 
чества просветляющих слоев, они остаются на 
уровне КнТ структуры 1, а затем ухудшаются. 

Поэтому мы предприняли попытку разделить 
20 нм слой серебра пополам (г/омй/е-структура). 
Расчёты показали, что это приводит к совершенно 
новому результату. При определенных комбинаци- 
ях геометрических и материальных параметров 
слоёв (см. структуры 2 и 3 на рис. 5) наблюдается 
заметное увеличение коэффициента пропускания 
и практически полное отсутствие отражения в об- 
ласти длин волн 450. ..750 нм (см. рис. 6). 

Оптимизация при дальнейшем увеличении ко- 
личества просветляющих слоев не даёт заметного 
преимущества перед структурами 2 и 3. 


Важным обстоятельством является то, что ко- 
эффициент отражения ТОП в ПК области практи- 
чески не меняется при переходе от структуры с од- 
ним слоем серебра (20 нм) к структуре двумя слоя- 
ми (10+10 нм) и составляет ~94 %. 


Ті0 2 37 нм 
Ад 20 нм 
ТЮг 35 нм 

стекло 4 мм 


8п0 2 48 нм 
Ад 10 нм 
Ті0 2 63 нм 
Ад 10 нм 
8п0 2 48 нм 
стекло 4 мм 


8п0 2 44 нм 
Ад 10 нм 
8п0 2 40 нм 
А1 2 0 3 9 нм 

8п 0 2 52 нм 
Ад 10 нм 
8 п0 2 45 нм 

стекло 4 мм 


Рис. 5 . Варианты структуры ТОП 


Несколько слов о технологических аспектах 
проблемы ТОП. Наши расчёты выполнены в пред- 
положении, что просветляющий слой выполнен из 
диоксида титана или системы 5п0 2 /А1 2 0 3 /8п0 2 . 


Объективно, по своим оптическим свойствам, Ті0 2 
способен обеспечить более высокое качество прос- 
ветления. Но скорость осаждения 3п0 2 в магнет- 
ронных диодах на постоянном токе в 5. ..10 раз вы- 
ше. Это обстоятельство может быть решающим 
при выборе технологической схемы. 

Необходимо помнить и о барьерных слоях, функ- 
циональное назначение которых - защитить серебря- 
ную плёнку от воздействия плазмы, содержащей ио- 
ны кислорода. Чтобы предотвратить влияние барьер- 
ных слоев на оптические свойства ТОП, они наносят- 
ся с минимальной толщиной. Поэтому их не учитыва- 
ли в процессе оптимизации. На первый взгляд, види- 
мо, достаточно наносить их с одной стороны слоя се- 
ребра (с той, которая подвергается воздействию ионов 
кислорода плазмы). Но имеет место фактор диффузии 
серебра в соседние слои, вызывающей деградацию 
плёнки. Его тоже надо принимать во внимание. 



Рис. 6. Коэффициенты пропускания и отражения различных 
структур ТОП (рис. 5) и чистого стекла 

Заключение 

1. Нами подтверждено, что эволюционный метод 
является эффективной процедурой оптимиза- 
ции оптических свойств теплоотражающих 
покрытий. 

2. Оптимизированы расчетным путем оптические 
характеристики многослойных теплоотражаю- 
щих покрытий типа "оксид - серебро - оксид" 
на стеклянной подложке при нормальном угле 
падения света. Установлено, что при разделе- 
нии плёнки серебра на два слоя с просветляю- 
щей оксидной прослойкой между ними ширина 
области высокой прозрачности в видимой части 
спектра увеличивается со 100 до 400 нм. 

3. Коэффициент отражения подобных покрытий в 
диапазоне длин волн 450.. .750 нм в сопоставлении 
с исходным стеклом снижается с 10 до 1 и менее %. 
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